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We establish that QED3 can possess a critical number of flavours, Nf, associated with dynamical 
chiral symmetry breaking if, and only if, the fermion wave function renormalisation and photon vac- 
uum polarisation are homogeneous functions at infrared momenta when the fermion mass function 
vanishes. The Ward identity entails that the fermion-photon vertex possesses the same property and 
ensures a simple relationship between the homogeneity degrees of each of these functions. Simple 
models for the photon vacuum polarisation and fermion-photon vertex are used to illustrate these 
observations. The existence and value of Nf are contingent upon the precise form of the vertex 
but any discussion of gauge dependence is moot. We introduce an order parameter for confine- 
ment. Chiral symmetry restoration and deconfinement are coincident owing to an abrupt change 
in the analytic properties of the fermion propagator when a nonzero scalar self-energy becomes 
insupportable. 
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I. INTRODUCTION 

Analogous to quenched quantum chromodynamics 
(QCD), quenched quantum electrodynamics in three di- 
mensions (2 spatial and 1 temporal - QED3) is confining 
because it has a nonzero string tension This feature 
persists in the unquenched theory if massive fermions cir- 
culate in the photon vacuum polarisation • That mass 
can either be explicit or dynamical in origin. The pos- 
sibility of dynamical mass generation in the chiral limit 
theory is of material interest because then QED3 has 
even greater qualitative similarity with QCD. 

Dynamical chiral symmetry breaking (DCSB) explains 
the origin of constituent-quark masses in QCD and un- 
derlies the success of chiral effective field theory [3j. 
Its importance in QCD has ensured that the possi- 
bility of DCSB in QED3 has received much atten- 
tion. The Dyson-Schwinger equations (DSEs) provide 
a natural framework within which to explore this and 
related phenomena. These methods are reviewed in 
Refs. 

Interest was stimulated by a finding that in QED3 
with Nf massless fermions a truncated system of DSEs 
yields a critical number of flavours, Np such that DCSB 
is only possible for Nf < Nf. A debate has subsequently 
ensued, which focuses on the reliability of the truncations 
employed and offers an alternative possibility; namely, 
that the mass function does not vanish at some criti- 
cal number of flavours but is instead exponentially sup- 
pressed with increasing Nf ■ Our analysis is intended 
to crystallise the issues and thereby eliminate the points 
of contention. 

It is worth remarking that lattice-regularised QED3 
can currently state only that N L f > 1, if it exists at all 



[111 ] . An impediment to reliable lattice calculations is the 
mass hierarchy feature of QED3; viz., any dynamically 
generated mass-scale is at least one order of magnitude 
smaller than the natural scale, which is set by the dimen- 
sioned coupling e 2 (see, e.g., Ref. [Hj]). 1 

In parallel with its relevance as a tool through which 
to develop insight into aspects of QCD, QED3 is also of 
interest in condensed matter physics as an effective field 
theory for high-temperature superconductors [HI, Q> Il5j 
and graphene [H,[l3]. Naturally, the nature of the theory 
at Nf — 2 is most interesting for true electronic systems. 

We describe some basic features of QED3 in Sec.HIl 
and introduce the coupled DSEs relevant to our analy- 
sis. Section IIIII canvasses the issue of gauge covariancc 
in truncations of the DSEs. It also elucidates neces- 
sary and sufficient conditions for the existence of a criti- 
cal number of flavours, in particular a collusion between 
Schwinger functions at infrared momenta. This is illus- 
trated using Ansdtze for the vacuum polarisation and 
fermion-photon vertex in Sec. lIVI In addition, we intro- 
duce a model-independent order parameter for deconfine- 
ment and demonstrate simultaneity of chiral symmetry 
restoration and deconfinement. We expect that this re- 
sult remains valid with all consistent truncations. Sec- 
tion [V] provides a summary and perspective. 



This is fortunately not the case in QCD, which possesses a di- 
mensionless but running coupling. Its evolution is determined by 
a mass-scale whose magnitude is characteristic of all dynamically 
generated mass-dimensioned quantities in the theory. 
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II. COUPLED EQUATIONS 

In three dimensions there are two inequivalent 2x2 rep- 
resentations of the Euclidean Clifford algebra {7^, "f u } — 
2<5 M „. Hence, two-component spinors are sufficient to de- 
scribe spinorial representations of the Lorentz Group. In 
this case, however, a mass term of any origin is odd under 
parity transformations. That can be avoided by employ- 
ing four-component spinors and a 4 x 4 representation of 
the Clifford algebra e.g., the set {71,72,74}, taken from 
four dimensional theories, with 75 := —717274. There 
are then two different mass terms that can be written in 
the Lagrangian; viz., "m/ipip" and "m-0i[7 3 , 75]^." The 
former is obviously analogous to the natural form in four 
dimensions and it is invariant under parity transforma- 
tions. We therefore define the theory through this term. 

In order to study confinement and DCSB we consider 
the gap equations for the fermion and photon; namely, 
in a theory with Nf fermions of mass m, 

Sip)- 1 = i 1 - P A(p) + B(p), (1) 
= [i 7 -p + M(p)]/Z(j>), (2) 
= 27 ■ p + m + £(p) , (3) 

S(p) = e 2 / D^(p-q) 7lM S(q)r v (q,p), (4) 
Jq 

where = j d?q/(2n) 3 and, with q± — q ± fc/2, 

D-l{k) = [^-(l-i/OM^ + n^fc), (5) 
n M „(fc) = [k 2 *„„ - k^K] n(fc) = : r^{k) k 2 n(fc),(6) 

= -N f e 2 f tv^S(q + )T v (q +> q_)S(q^).(7) 

Jq 

NB. The massless theory is straightforwardly defined by 
setting the Lagrangian mass m — and the quenched 
theory is obtained by writing II(fc) = 0. 

These equations are written in a general covariant 
gauge, fixed by the parameter £. Plainly, they are cou- 
pled to each other and to equations that we have not 
written explicitly; e.g., that for the fermion- photon ver- 
tex, r„. In order to draw reliable conclusions from the 
analysis of these equations a systematic and nonpertur- 
bative truncation scheme should be employed. One such 
scheme was introduced in Refs. (l8l. [l9j. 

The coupling a = e 2 has mass-dimension one. Fur- 
thermore, QED3 is super-renormalisable and therefore, 
in contrast to QCD, no ultraviolet divergence can arise 
whose regularisation would introduce a new mass-scale. 
Hence, a sets the scale in the massless theory and there 
is a direct connection between a and the scales of con- 
finement and DCSB. 

With the parametric dependence of the Schwinger 



functions displayed explicitly, these features entail 

77 77? 

aS(p;m;a;N f ) = S(^; -; l;N f ), (8) 
a a 

D 777 

n(p; m; a; Nf) = n( — ; — ; 1; Nf), (9) 
a a 

T ll {p,q;m;e 2 ;Nf) = r M ( — , — ; — ; 1; Nf); (10) 

a a a 

i.e., that the pointwise evolution of the theory's 
Schwinger functions at an arbitrary value of a can be ob- 
tained through scaling of their a = 1 behaviour. In con- 
sequence there cannot be a phase transition associated 
with changes in a: whether QED3 is confining and/or 
exhibits DCSB is independent of the coupling's value. 
On the other hand, the state of the theory can respond 
to changes in the dimensionless parameter Nf. 

III. DCSB - THEORY 

In bald terms DCSB is the existence in the chiral 
limit of a M(p) ^ solution to Eq. Q. Such a solu- 
tion expresses the generation of mass from nothing solely 
through interactions. 2 In order to determine whether 
that is possible one can begin by analysing the gap 
equations at leading order in the truncation scheme of 
Refs. [H, which corresponds to using 

r M fr,<?) = 7„- (11) 

In this "rainbow-truncation" Eqs. and (JT]) form a 
closed set, which has long been the subject of scrutiny; 
see, e.g., Refs. 0,113 and references therein. 

A. Landau gauge 

Before reporting the results of such an analysis, we 
consider it important to discuss whether there is a spe- 
cial choice of gauge parameter for which Eq. (fTTj) is most 
sound. If so, then one can couple that particular gauge 
choice with Eq. (JTTJ) and define the vertex in any other 
gauge via the Landau-Khalatnikov-Fradkin (LKF) trans- 
form of Eq. pT|) [2ll . |22| . |23j . This procedure was advo- 
cated in Ref. Q because it guarantees gauge covariance. 

Landau gauge, £ = 0, occupies this special place. In 
Landau gauge the one- loop contribution to A(p) vanishes 
in any number of dimensions [24j |. Owing to the sim- 
plest Ward-Takahashi identity, this feature amplifies the 
domain whereupon Eq. (TTTj) can be described as a point- 
wise accurate approximation. Additionally, in the con- 
text of the reasonable vertex Ansdtze [25( that have been 



2 Whether the vacuum associated with this solution is preferred 
over that corresponding to the M(p 2 ) = solution; viz., has the 
higher pressure, can subsequently be determined. It certainly 
docs in QCD [2(J. 
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employed in DSE studies, this means that sensitivity to 
model-dependent differences between the forms is least 
noticeable in Landau gauge. Moreover, whenever it is rel- 
evant, £ = is a fixed point of the renormalisation group; 
i.e., the gauge parameter is momentum- independent. In- 
deed, it is always zero. 

Numerous studies have tried to construct a model ver- 
tex that guarantees gauge covariance of the fermionpror. 
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1361 . 137| and the photon vacuum polarisation @, 
QED3 and QED4. The lesson in the observations above 
is that such attempts are misguided. Suppose a Landau 
gauge Ansatz is truly correct, then its form in any other 
gauge is precisely reproduced by the LKF-transform. If, 
on the other hand, the Ansatz is wrong, then its in-built 
dependence on the gauge parameter will not be better 
than that which is obtained by defining the vertex in 
different gauges via the LKF transform. Hence there is 
nothing to be learnt from exploring the gauge parame- 
ter dependence of quantities obtained with a given vertex 
Ansatz. 3 Our view is that irrespective of the choice made 
for r^, and there are good reasons to employ something 
more sophisticated than Eq. (fTTj) . the Landau gauge re- 
sults define the model entirely. 



B. Conditions for a critical number of flavours 

In Ref. 0], in conjunction with free massless fermion 
propagators, S(p) = — i/7 • p, Eq. (fTT|) was used in 
Eqs. © and ® to obtain 



IL(k)=N f a- 



(12) 



This result was subsequently inserted into Eq. ([4]) via 
Eq. iJSJ) and, subject to the further truncation A(p) = 1, 
an equation for M(p) obtained; viz., with a = 1, 



M(p) 



2ir 2 p 



dq 



qM(q) 

q 2 + M{ q y 



In 



p + q + 



\p-q\ + ^ 



(13) 

This procedure was described as providing the mass func- 
tion at leading-order in a 1 /^/-expansion and yields the 
result that DCSB is only possible on the domain 



N f < Nf 
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3.24. 



(14) 



This claim opened a debate, described in Sec. 3 of Ref. [4(. 
The other widely canvassed possibility is that the mass 



function does not vanish but is instead merely exponen- 
tially suppressed with increasing Nf [l(| • One aim of the 
following analysis is to eliminate the points of contention. 

Given that H(k) must be a monotonically decreasing 
function in QED3, our first observation is that a critical 
number of flavours, Nf, can only exist in the massless 
theory if in the neighbourhood of this value 



1 



n(fc) 



N f ~Nf 



0. 



(15) 



To establish this, consider the opposite; namely, 
U(Q;Nf ~ Nf) = a, for some a € (0, 00). Under the 
conditions just described 



1 > 



1 



> 



1 



1 + n(fc) ~ 1 



> , Vfc > 0. 



(16) 



Equation ([4]) can only support a gap if the kernel pos- 
sesses sufficient support. A lower bound on the support 
is obtained by employing the replacement: II(fc) — > a. 
The equation thus obtained is actually the same as the 
quenched equation except for a rescaling of the coupling; 
viz., a — > a — cr/(l + a). In the quenched theory 
the coupling serves merely to set a mass-scale and has 
no effect on whether or not chiral symmetry is dynam- 
ically broken. Thus, if chiral symmetry is not broken 
for 11(0; Nf ~ Nf) = a, then there exists a N' f < N c f 
such that 11(0; Nf) = a 7^ a and chiral symmetry is still 
unbroken. Therefore, Nf cannot be a critical number 
of flavours. Hence no monotonically decreasing function 
satisfying Tl(k) < a < 00, Vfc > can generate a critical 
number of flavours. 

If a theory possesses a critical number of flavours, then 
in the chiral limit and for Nf ~ Nf one can neglect in the 
Schwinger functions all instances of feedback associated 
with B(p) ^ 0. Hence Eqs. (j4]) and ([7]) can be studied in 
the forms 



B(p) = f D^(p-q) 
J a 



J^ltr^Tyfap), (17) 



P 2 A(p) = P 2 + 



l D ^ p - q) 7W) 



1 



x-tr (-i-y ■ phui-h ■ Q)T v (q,p), (18) 



n(fc) 



N 



1 



xtr [7 M (-«7 • q+)r v (q+, q-)(-i-f ■ ?-)] .(19) 



2fc 2 J q qiA(q + )q±A(q-) 



3 In this connection it is materially important to bear in mind that 
no vertex Ansatz can be consistent with perturbation theory 
if its expression in other than Landau gauge omits an explicit 
dependence on the gauge parameter; i.e., the correct vertex must 
possess a dependence on £ ^ in addition to that implicit in 
the fermion propagator which arises through resolving the Ward- 
Takahashi identity |33|1 . 



(NB. Owing to Eqs. (J8]) - (fT0|) . here and hereafter we set 
a = e 2 = 1.) In deriving Eqs. (fTTJ) - (TT9"1) we use d the 
fact that without feedback; i.e., for B(p) = 0, all terms 
in r M contain a single power of the Dirac matrices. No 
term is present that can drive the appearance of another 
structure. This entails that the equations for A(p) and 
Tl(k) form a closed, coupled pair, one whose solutions will 
only ultimately influence the result for B(p). 
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Consider Eqs. ([18]) and (|T9|) . Equation (fl2|) describes 
the polarisation for p » 1, on which domain one also 
has A(p) -1«0. However, it is the behaviour of these 
functions for infrared momenta, p -C 1, that is crucial 
to the existence or otherwise of a gap. Reference [25[ 
indicates the possibility that in the absence of an explicit 
or dynamically generated fermion mass, A(p) can possess 
homogeneous power-law behaviour in the infrared whilst 
nevertheless satisfying its perturbative ultraviolet limit. 
We therefore consider the implications of 



A(p) = a p S , S > 0, p< 1. 



(20) 



Homogeneity entails A((p) = C 8 A{p) . From this, the 
absence of a fermion mass and, importantly, the Ward 
identity, it follows that 



(21) 



rp(Cp,C«) = C*r p (?,<?), P ,q<i. 



Subject to these considerations, inspection of Eq. (|19|) re- 
veals that infrared momenta dominate the integration do- 
main therein. Therefore 



n(Cfe) = c 



II(k) , k < 1 . 



(22) 



NB. Dominance of infrared momenta cannot be argued 
for S < 0. Moreover, Eq. ([22]) is consistent with Eq. (fT5|) . 

Equation (fT8|) has a constant term on the left-hand- 
side. This is plainly the driving term associated with the 
ultraviolet behaviour A(p) = 1, p ^> 1, and therefore can- 
not have any impact on the infrared behaviour of A(p), 
a result illustrated clearly in Ref. [37|- Hence, focusing 
on momenta within the infrared domain, we substitute 
Eqs. ([20 )1 - (p2] ) into Eq. (18]) and find 



A((p)=( s A(p); 



(23) 



namely an internally consistent result. The value of 5 
depends on the exact form of the dressed- fermion-photon 
vertex [33], which is not completely constrained by the 
Ward-Takahashi identity. 

It is thus evident that if there exists a critical number 
of flavours, Nj, then for Nj ~ Nj the fermion wave func- 
tion renormalisation function and the photon vacuum po- 
larisation are, for infrared momenta, homogeneous func- 
tions. They have degree S and [— (1 + 8)], respectively, a 
particular relationship that is enforced by the Ward iden- 
tity. The converse can also be argued from Eqs. (0| and 
([71) . If there exists a value of Nj such that for Nf ~ Nj 
A(p) and n(fc) are homogeneous functions at infrared mo- 
menta, then at this value of Nj chiral symmetry is not 
broken. This is plain because DCSB would introduce an 
infrared mass-scale, /i, and the wave function renormal- 
isation and vacuum polarisation would saturate at some 
constant value for momenta less than this scale. 

It follows that there can exist a critical number of 
flavours for DCSB iff the fermion wave function renor- 
malisation function and the photon vacuum polarisation 



are homogeneous functions at infrared momenta. 4 

These results have an important impact on Eq. {T7J); 
namely, for any value of 8, 



B(Cp) = B(p) ,p<l 



(24) 



An infrared collusion is apparent in the theory: the in- 
frared anomalous dimensions of the wave function renor- 
malisation and vacuum polarisation cancel from the 
equation for the scalar self energy and the scaling be- 
haviour of the integrand is £°. Once again, the Ward 
identity driven result in Eq. ([2"Tj) is essential for this to 
occur. Equation (|24[) states that it is consistent with 
homogeneity, Eqs. (|2"Tj) - ([2"3"| . for the scalar part of the 
fermion self energy to assume a constant value - possibly 
nonzero - at infrared momenta in the neighbourhood of 
a critical number of flavours, Nf. 

To underline the importance of this outcome, sup- 
pose that the fermion wave function renormalisation, 
fermion-photon vertex and photon vacuum polarisation 
are treated inconsistently; e.g., such that in Eq. ([22)) the 
exponent is [—(1 + 6')]. Then the integrand in Eq. (JT7J 
scales as C, s ~ s , in which case the only admissible solution 
is B(p) = 0. Such a result precludes the existence of a 
critical number of flavours. This analysis reiterates and 
elucidates the observation of Ref. [3l[ . 



IV. ILLUSTRATION 
A. DCSB 

In order to provide a straightforward illustration via 
the fermion gap equation, Eq. (g]), of the results discussed 
in Sec. lIIIl we adapt the vacuum polarisation calculated 
from numerical solutions of the gap equation in Ref. 0] ; 
namely, 



U{p,q;£ v )=V{p,q;£ P )U{ P -q), 



H(fc) = Nf 



1 



8 



+ <;bc 



—ck 



(25) 
(26) 



In Eq. ([H]) we use a = 0.20, b = 0.088 and c = 7.8, 
which owe their nonzero values to DCSB and were cal- 
culated for Nf = 1. NB. For a = = b, Eq. (2BJ) 
reduces to Eq. (fT2"[) : i.e., the leading order result in a 
1 /A^j-expansion. 

In adapting this form to our purposes we have intro- 
duced two new factors. The quantity 



-2(N f -l)/p 



M(p = 0; Nf) 
M(p = 0;N f = 1) 



(27) 



4 In the fermion gap equation one can represent a Higgs mass, r/, 
for the photon via an infrared singular polarisation: Iljy(fc) = 
r] 2 /k 2 . Hence the discussion above also applies in that case. 
Reference '381 can be viewed as a numerical verification of this. 
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allows a suppression of M(p) with increasing Nf to feed 
back into the vacuum polarisation. This enables a sce- 
nario in which Eq. (|15|) can arise dynamically in concert 
with the appearance of a critical number of flavours. We 
have also included the factor 



V{ P ,q;£ v ) = 



1 



1 



A(p) A(q) 



(28) 



the presence of which permits the model to realise 
Eq. (|22[) dynamically when l-p = 1. Similar reasoning 
leads us to write 



(29) 



with 



v^M = {\ [A{p) + A{q)h % = = l • 00) 

At this point our illustrative model for the gap equation 
is completely specified. 

We solve Eq. Q with a — e 2 = 1 , the Landau gauge 
photon propagator obtained from Eqs. ([5]) and ([26]) - ((281) 
and the fermion-photon vertex in Eqs. ([2"9"| . ([50)1 . The up- 
per panel of Fig.[T] describes the nature of DCSB in the 
model. It is apparent that a critical number of fermion 
flavours exists if, and only if, the vacuum polarisation 
and vertex dressing conspire at infrared momenta to com- 
pletely eliminate the influence of the fermion's vector self- 
energy from the equation for the scalar self-energy. Oth- 
erwise, M(p; Nf) is merely suppressed exponentially with 
increasing Nf. 

The solid curves in Fig.Q] exhibit a sudden drop at 
Nt = 1.55. This can be understood via Fig. [2J which 
depicts the vacuum polarisation of Eq. (f2"o]l . As we ad- 
vertised, the polarisation acquires a dependence on Nf 
because a suppression of M(p) feeds back into the po- 
larisation through This mocks up a simultaneous self- 
consistent solution of the fermion and photon gap equa- 
tions. The domain on which H(k) differs from the 1/iV/- 
result in Eq. (| 1 2j> shrinks with increasing Nf, until at 
Nf = 1.55, owing to the calculated self-consistent feed- 
back, it no longer contributes measurably to the fermion 
gap equation. Thereafter, this gap equation's solution 
also follows the l/7V/-result. 

The following observations are relevant. It is II(p, q; 1) 
in Eq. ([25]) that self-consistcntly realises Eq. (|22|) . Also, 
in the upper panel of Fig.[U which depicts p, the solid 
and short-dash-dot curves do not lie upon one another 
because of the normalisation factor in Eq. ((2"T)) : it is plain 
from the figure that our Ansatze yield a Nf = 1 result 
for the mass function which is larger than that obtained 
in a 1/Nf expansion. The curves are indistinguishable if 
for Nf > 1.55 we normalise our solution via the Nf = 1 
result obtained using the 1/Nf expansion. 

It is now clear that in the neighbourhood of chiral sym- 
metry restoration Eq. ([2"o]) is identical to Eq. (fT2")l , and 
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FIG. 1: Upper panel - Calculated evolution of the dimension- 
less ratio p in Eq. 1)270 with the number of flavours, Nf. dotted 
curve - (lv =0,iv = 0) in Eqs. (f26|) . J29J, respectively, which 
represents the complete suppression of feedback into the gap 
equation from the solutions of other DSEs; short-dashed curve 
- (lv = l,lv = 0); long-dashed curve - (lv =0,£v — 1); an d 
solid curve - (lv = l,^v = 1)- Only the last case, which al- 
lows a dynamical modification of the photon vacuum polarisa- 
tion and fermion-photon vertex to influence the gap equation, 
exhibits a critical number of flavours: Nj — 3.24, Eq. (|31)1 . 
Short- dash- dot curve - For comparison, the result obtained 
via a l/iV/-expansion. Lower panel - The confinement order 
parameter n(Nf), Eq. (|35|l . Simultaneity of chiral symmetry 
restoration and deconfinement is a robust result. 



hence our Ansatze yield 



Ni 



32 

72 



3.24 and S(N c f ) 



1 

12 



0.083. (31) 



We re-emphasise that these values are typical but not 
definitive. Extant studies suggest that a more sophisti- 
cated treatment of the complex of relevant, coupled DSEs 
would yield results that differ by < 25%. 
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FIG. 2: The polarisation in Eq. ([26} plotted as k Tl(k) and 
calculated for a range of values of Nf. dot-dashed, Nf = 1.2; 
short-dashed, Nf = 1.4; solid, Nf = 1.5; and dotted, Nf > 
1.55. The ^-dependence owes to ? in Eq. (|27[> and is self- 
consistently determined. Even with the much magnified axes 
scales in this plot, for Nf > 1.55 the results are indistin- 
guishable from each other: they lie upon the 1/Nf -result in 
Eq. HU). 



B. Confinement 

Quenched QED3 is confining because it has a nonzero 
string tension [ij . This feature persists in the unquenched 
theory but only so long as 11(0) is finite [2j. As we have 
seen, 11(0) is finite iff the fermions are massive: in QED3 
only massless fermions can screen completely. 

The implications of this for the fermion 2-point func- 
tion can be read from Sec. 2 of Ref. If one writes 



S(p) = -17 • po- v (p 2 ) + a s (p 2 ), 
then in the absence of confinement (x = p 2 ) 
d 2 



dx 2 



crv(x) > , Vx > . 



(32) 



(33) 



On the other hand, S(p) describes a confined excitation 

d 2 



if 3x n > : a v(x) 



0. 



(34) 



These statements are associated with the realisation of 
confinement through a violation of the axiom of reflec- 
tion positivity. Any 2-point Schwinger function with an 
inflexion point at x = p 2 > 0, Eq. (|34[) . must breach the 
axiom of reflection positivity. This entails that the asso- 
ciated elementary excitation cannot appear in the Hilbcrt 
space of observables. 

In Fig. [3] we depict the evolution of a' v (x) calculated 
in the model described in Sec. IIVI with three values of Nf 
below Nf and another above. For Nf < N*, chiral sym- 
metry is dynamically broken. The conspicuous maximum 
in (— a'y) for Nf = 1, 1.5, 2 signals fermion confinement in 
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FIG. 3: Calculated evolution of ln(—a' v ): solid curve - cal- 
culated with Nf = 1; short-dashed curve - Nf = 1.5; long- 
dashed curve ~ Nf — 2; and dashed curve - Nf — 4. The 
conspicuous maximum in ln(— a'y) for Nf = 1, 1.5, 2 signals 
fermion confinement in the domain of Nf whereupon chiral 
symmetry is dynamically broken. (Recall that e 2 = 1.) 



the domain of Nf whereupon chiral symmetry is dynam- 
ically broken. The location of the minimum in o~' v (x), 



migrates to x = as Nf 



N 



f 



One can iden- 



tify x c as an order parameter for deconfinement. NB. 
The location of the minimum positioned farthest from 
x = p 2 = should be used. In an asymptotically free 
theory, if there are any minima at all, then there will be 
one most distant. 

It is evident in the lower panel of Fig. [T] that the quan- 
tity 



•XNf) 



XNf = 1) 



(35) 



exhibits precisely the same behaviour as p(Nf). The fig- 
ure demonstrates that when internally consistent feed- 
back is enabled, chiral symmetry is not dynamically bro- 
ken and the fermions are not confined for Nf > Nf. 
Plainly, the chiral symmetry restoration and deconfine- 
ment transitions are simultaneous in our illustrative 
model. The basic causal connection is a dramatic change 
in the analytic properties of the propagator which ac- 
companies the disappearance of a nonzero fermion scalar 
self-energy. 

It is noteworthy that DSE studies of QCD, in rainbow- 
ladder truncation with the model kernel of Ref. [39(1 , find 
the same results for these transitions at T ^ [4fj] and 
(j, ^ Q [41]. In addition, Ref. [12] describes a model that 
exhibits a line of simultaneous transitions in the physi- 
cal quadrant of the (T, /x)-plane. Indeed, deconfinement 
and chiral symmetry restoration are coincident in all self- 
consistent studies of concrete models of continuum QCD 
that exhibit both phenomena, and in numerical simula- 
tions of lattice-regularised QCD [43j. The methods de- 
scribed herein could be employed to extend the study of 
QED3 at nonzero- (T, /i) in Ref. (44[, and explore whether 
and under which conditions chiral symmetry restoration 
and deconfinement are coincident. 
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V. APERQU 

We argued that QED3 can possess a critical number 
of flavours, Nf, above which dynamical chiral symme- 
try breaking is impossible if, and only if, in the neigh- 
bourhood of Nj the fermion wave function renormalisa- 
tion and the photon vacuum polarisation are homoge- 
neous functions at infrared momenta. The Ward iden- 
tity guarantees that the fermion-photon vertex also pos- 
sesses the homogeneity property and ensures a simple 
relationship between the homogeneity degrees (anoma- 
lous dimensions) of each of these functions. One cannot, 
however, conclude from this analysis that QED3 must 
possess a critical number of flavours. The existence and 
value of Nj are contingent upon the precise form of the 
fermion-photon vertex. Notwithstanding this, it is note- 
worthy that all extant studies which employ a reasonable 
vertex Ansatz and treat QED3's gap equations consis- 
tently yield a value Nf ~ 3.5 . It is therefore almost 
certain that chiral symmetry is dynamically broken for 
Nf = 2 QED3, the number relevant in condensed matter 
applications of the theory. 

While our analysis does not establish the existence of 
Nf it does obviate any further debate on the question 
of whether the fermion mass vanishes at some finite Nf 
or is merely exponentially damped with increasing Nf. 
We have elucidated the conditions under which a criti- 
cal number of flavours can exist and the mechanism by 
which it then arises; viz., infrared collusion. All extant 
studies that report an exponential suppression of M(0) 
have suppressed that mechanism. 

Our discussion should also forestall any further con- 
sideration of the gauge dependence of Nj and other 
physical quantities. We have made plain that Landau 
gauge occupies a special place in gauge theories. It is 
the gauge in which any truncation or sound Ansatz for 
the fermion-photon vertex can most reasonably be de- 
scribed as providing a pointwise accurate approximation. 
The vertex in any other gauge should then be defined as 



the Landau-Khalatnikov-Fradkin transform of the Lan- 
dau gauge Ansatz. The sensible implementation of this 
procedure guarantees gauge covariance and hence renders 
moot any question about the gauge dependence of gauge 
invariant quantities. 

In order to illustrate our arguments we employed a 
simple model for the photon vacuum polarisation and 
fermion-photon vertex. Within this model we calculated 
Nf and the infrared anomalous dimension. Furthermore, 
we introduced a new, model-independent order parame- 
ter for the confinement of elementary excitations. In our 
model deconfinement and chiral symmetry restoration 
are coincident at Nf. This owes to an abrupt change in 
the analytic properties of the 2-point fermion Schwinger 
function when a nonzero scalar self-energy becomes im- 
possible. It is plausible that this mechanism and result 
persist in QED3 proper and, indeed, in general. 

We remark in closing that all the lessons learnt in this 
study are relevant to QCD. The simplest points are that 
Landau gauge is of particular utility; and confinement 
and DCSB might in general be tied to, and connected via, 
the analytic structure of the fermion propagator, which 
itself is strongly influenced by the possibility of a dynam- 
ically generated scalar self-energy. The phenomenon we 
have labelled as infrared collusion could be relevant to the 
phase diagram of QCD in the plane formed by the cou- 
pling and the number of quark flavours; e.g., Ref. [lol |46T| . 
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